Mass correction and gravitational energy radiation in 
black hole perturbation theory 

Reinaldo J. Gleiser 



Facultad de Matematica, Astronomi'a y Ffsica, Universidad Nacional de Cordoba, 
Ciudad Universitaria, 5000 Cordoba, Argentina. 

Abstract 

Using second order black hole perturbation theory, we show that the difference 
between the ADM mass and the final black hole mass, computed to the lowest 
significant order, is equal, to the same order, to the total gravitational radiation 
energy, obtained applying the bandau and bifschitz (pseudotensor) equation to 
the first order perturbation. This result may be considered as a consistency check 
for the theory. 



PACS numbers: 04.25.Nx, 04.70.Bw 



1 Introduction 

The analysis of the motion of gravitating bodies in the strong field region is, clearly, 
one of the most important problems in General Relativity, since it is here where one 
would expect to find the most pronounced effects, allowing, in principle, for a clear 
confrontation of the theory with observation. Perhaps the most important of these 
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effects, and the one to which the largest theoretical and experimental efforts is being 
devoted nowadays, is the emission of gravitational radiation through gravitational waves. 
As is well known, however, and in contrast with the situation in other theories, e.g. 
electromagnetism, the whole subject of gravitational waves and gravitational radiation in 
the context of general relativity is a very complex one, both from the point of view of their 
physical interpretation as from that of the difficulties in constructing solutions, (either 
exact or sufficiently approximate), appropriate for the modeling of relevant astrophysical 
systems. 

From very general arguments, one expects that the most prominent sources of bursts 
of gravitational radiation, are the coalescences of two astrophysical objects, leading to 
a single black hole as the final result. These same general ideas indicate that most of 
the radiation, and therefore of the information to be observed at a large distance from 
the region of emission, should come from the last stages of this coalescence, where the 
system is close to, and rapidly approaching the final black hole stage. This has lead, 
among other lines of attack, to the idea that, given adequate initial data, one could 
consider these last stages as corresponding to the evolution of the perturbations of a 
suitable chosen black hole. 

A successful application of this idea, attributed to Smarr, was carried out by Price 
and Pullin , in the case of the head-on collision of two black holes, in the close approx- 
imation, applying a technique, originally developed by Regge and Wheeler 0,and Zerilli 
0, (see also ||), devised precisely for the analysis of linearized perturbations of black 
hole. By successful we mean here that the results given perturbatively in JIJ, are essen- 
tially identical to those obtained by "exact" numerical solution of Einstein's equations, 
for the same initial data, for a suitably chosen range of the parameters characterizing 
the problem. 

It should be clear, however, that if the numerical results had not been available, it 
would have been essentially impossible to asses the range of validity of the perturbative 
results, not only because it relays only on linearized theory, but also because of the 
intuitively appealing, but non rigorous nature of the identification of what is meant by 
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the "gravitational energy" radiated by the system. 

An important step in the direction of defining a range of validity for the perturbative 
treatment was given in ||, where the formalism was extended to second order, with the 
idea that second order corrections should serve as a sort of "error bars" on the first order 
computations. This formalism was applied to the head-on collision of black holes in ||, 
showing a very remarkable agreement between the "exact" result, the first order results 
and the associated "error bars". 

A somewhat different application of black hole perturbation theory, taken to second 
order, would be to analyze, instead of the gravitational radiation, the change in the mass 
parameter of the black hole. From a physical point of view, we expect that if an isolated 
astrophysical system contains, at a given time, certain mass-energy, and that if, after a 
long time, part of this mass-energy is radiated to infinity, the sum of the mass-energy 
remaining in the system and that radiated to infinity, should be equal to initial mass. 
In more rigorous terms, the initial mass-energy, for an isolated system, is given by the 
ADM mass. Since the process which we are envisioning is one where the final result is 
a black hole, the final mass-energy is simply that of the final black hole. The energy 
radiated to infinity should then be equal to the difference between these two masses. 
This comparison would provide an independent check of the validity of the computation 
of the radiated energy, as given in, e.g., [T|]. 

In this note we consider the problem of computing the correction to the mass param- 
eter in the lowest non trivial order in black hole perturbation theory. After reviewing 
briefly the formalism, we derive a general equation for the relevant metric perturbation 
function, and, under some general restrictions, and using the Zerilli equation, obtain 
the expression for the mass correction, in terms of the Zerilli function. It is remarkable 
that the resulting expression is identical to that obtained from the Landau-Lifschitz 
pseudotensor for the total gravitational energy radiation, as given in although this is 
not mentioned anywhere in the derivations. In some sense, we might even consider this 
reassuring proof of consistency of the perturbation expansion, as a "derivation" of the 
radiation formula, since, as will be seen below, only general properties of the solutions 
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of Einstein's equations are used in the proofs, without any reference to asymptotically 
flat or radiation gauges. We should, however, make clear that it is not the intention of 
this article to delve into the difficult problem of the definition of the energy in general 
relativity, and of the relative merits of the different "pseudotensors" or "complexes" that 
have been discussed in the literature |§ . In particular, the comparison is made only to 
the expression for the integrated power, and it is known that different prescriptions may 
give here the same final answer. We also remark that the equivalence between the Bondi 
energy flux and that given by the Landau-Lifschitz complex, in suitable asymptotically 



flat space-times, has been shown in [10 



2 The perturbative expansion 

The essence of the perturbation method is the expansion of the metric in the form 

9 l u,(e)=g® + eg$+<?gV> + 0(<?) (1) 

where g^e) represents a one parameter family of solutions, and gj® some known exact 
solution of Einstein's equations, leading to an expansion of the Einstein equation into 
a hierarchy of equations, ordered also by the parameter e, where the solution of each 
order, considered as an initial value problem, requires the solution of all previous ones 
0. Some relevant issues, including gauge invariance, are discussed in more detail in 0. 
(See also §.) 

In our case, g@J corresponds to the Schwarzschild black hole, and e to some parameter 
characterizing the departure of the initial data from that of a black hole. A particularly 
suitable framework for this problem is given by the Regge- Wheeler M formalism, where 
the perturbations to any given order are given as a multipolar expansion in the angular 
variables 6, and (p. The multipoles are, in turn, separated into even and odd type, for any 
given L, the order of the multipole. The L = 0, (monopole), terms contain information 
on the mass of the system. We shall assume that the leading (first order in e) part of the 
perturbation corresponds to the L = 2, even, (quadrupole) terms. (This is, for instance, 
the case considered in Therefore, all other terms are, at least, of order e 2 
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We choose from the outset a Regge - Wheeler gauge 0. This means that the non 
vanishing metric metric coefficients are 

9u = -(l-^)[l-eH^\t,r)P 2 (9)-e 2 H^P 

9rr = (l- 2 ^)~\l + eH?{t : r)P 2 {9)+e 2 H^P Q 
g ee = r 2 [l + eK^(t, r)P 2 {6)} (2) 
g H = r 2 sm(9) 2 [l + eK^(t,r)P 2 (9)] 

where the Pl, L = 0, 2 are Legendre polynomial in cos(#). The factor e makes explicit 
the perturbation order. The functions #g (t, r), Hi(t,r), H^it^r), and K^(t,r) 
characterize the first order, L = 2, even, perturbations, while, Hq (t,r), and H^\t,r) 
describe the L = 0, second order perturbations. The general second order expansion 
contains terms also for L = 2 and L = 4, but these are decoupled, to this order, from 
each other, and from the L = terms, and we shall not consider them here. 

We remark that although, in general, the L = perturbations would also contain 
terms of the form 

gtr = e 2 H[ 2) (t,r)P (9) 

gee = r 2 e 2 K^(t,r)P (9) (3) 
g H = r 2 sm(6) 2 e 2 K^(t,r)P (e) 

it is always possible to choose a gauge where H[ 2 \t,r) = 0, and K^ 2 \t,r) = 0, so 
we make that simplifying choice. Furthermore, with this choice, if we assume that the 

(2) 

first order perturbations vanish, the solution of Einstein's equations for if^ is time 



independent, and of the form 



r — 2M 



where C is a constant. 

Now, if we consider as the only "perturbation" a change in mass 8M, it can be seen 
that 

rr(2) 25M 



Therefore, we identify the constant C with twice the correction to the mass. 

Taking into account the different powers in e, and in a manner analogous to that in 
||, the Einstein equations for (§) separate into a linear set of equations for H^\t,r), 
Hi(t, r), H%(t, r), and K^\t, r), corresponding to terms linear in e, and a set of equa- 
tions linear in HQ 2 \t,r), and H^it^r), but containing "source"-like terms, quadratic in 
the first order functions. 

Regarding the first order functions, we recall that the general solution to the first 
order equations can be given in terms of the Zerilli function. Namely, any set ifg 1 ' ) (t,r), 
Hi(t,r), H^ftjr), K^(t,r) of solutions of these equations can be written in the form 



if (1) (t,r) 
HP(t,r) 

H^(t,r) 
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r 2 + rM + M 2 
r 2 (2r + 3M) 



r J 



d_ 

dr 
2r 2 



2r 2 - 6rM - 3M 2 
r(2r + 3M) 



?/> (1) (t,r) + (r - 2M) 



dr 



dr 



6rM-3M 2 dip^(t,r) d 2 ip^(t,r) 



(r-2M)(2r + 3M) 
H?\t,r) 



dt 



drdt 



where ip^(t,r) is a solution of the Zerilli equation: 

d 2 ^(t,r) d 2 ij^\t,r" 



dr 



*2 



-K(r*)^ (1) (t,r) = 



dt 2 



where 



and 



V{r) 



r + 2Mln[r/(2M) - 1] 



M\ 4r 3 + 4r 2 M + 6rM 2 + 3M 3 
6(1-2— — 
r / 



(6) 

K^(t,r)(7) 

(8) 
(9) 

(10) 
(11) 



(12) 



r 3 (2r + 3M) 2 

We remark that, for sufficiently large r, ( |T0"D approaches the wave equation, and, 
therefore, for sufficiently small initial data, ip^(t, r) approaches a function of t — r*, for 
large r and t. 

Going now to the second order in e, L = 0, perturbations, we find that there are four 
non trivial equations (written as R^ v = 0, where R^ v is the Ricci tensor), for the second 
order functions. The equation R tr = takes the form 

dH?\t,r) 



dt 



s 



(13) 



6 



where S is quadratic in the first order perturbations. We assume a choice of initial data, 
(for t = 0), such that M is equal to the ADM mass, and, therefore, we should have 



lim 

r— >oo 



r^ 2) (0,r) 







(14) 



so that there is no correction to the mass for t = 0. 

Since on the other hand, for large t, the solution should approach the static black 
hole configuration, we should have 

25M 



and we find that 



and, then, 



lim 

t— >oo 



9H?\ t , r ) J6M _ m 
o dt r-2M 2 v ' ; 



Hi 2 \t,r) 



r -2M 



8M =- lim 

2 r^oo 



r / Sdt 

o 



(15) 



(16) 



;i7) 



Using the asymptotic properties of the Zerilli function, and a fair amount of algebra, 
(some details are given in the Appendix), we find, 
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2 


lim 






dt 


r— »oo 


To 


Jo 


dt 



which expresses the change in mass in terms of the first order perturbations, through the 
Zerilli function. But the right hand side of this equation is also immediately recognized 
as precisely (minus) the total gravitational radiation energy, computed from the Landau 
- Lifschitz (pseudotensor) equation [|]. 



3 Comments 

The main result of the present analysis, indicated by Eq. ([18]) is interesting in several re- 
spects. To begin with, we remark that it was obtained in the framework of second order 
perturbation theory, namely, the presence of the "source terms", containing the contri- 
butions from the lower (first) order perturbations was crucial in its derivation. Second, 
we find an expression for the radiated energy, working only in the Regge- Wheeler gauge, 
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that is totally independent of any consideration of asymptotically flat or radiation gauges 
and, or identifications of gravitational wave amplitudes. Furthermore, the expression, 
although computed in second order perturbation theory, contains only information from 
the first order quantities. Thus, in a sense, it may be considered as a "derivation" of 
the radiation formula to be used in first order perturbation theory. That this equa- 
tion coincides with the Landau-Lifschitz pseudotensor prescription, may then be taken 
as a reassuring (although partial) proof of the physical consistency of the black hole 
perturbation treatment. 

Finally we remark that the proof in this note has been limited to a restricted, al- 
though relevant, set of perturbations, since only the first order, L = 2 even case was 
considered. This choice was made mainly for simplicity and we expect that an analogous 
result should hold in the general case. This will be considered elsewhere. 



Appendix 

We include in this Appendix some computational details. First we notice that, from 
R tr = 0, after some simplifications using the first order equations, we find 

rH?\, = t - \ (j^,+sP*f\,-aP«f\,) 

J. K m H m - r -^L H ?H? + f -HP„ /f(") . (19) 

Assuming an outgoing wave boundary condition for large r, (appropriate for pertur- 
bative initial data), one can show that ip(t, r) admits the following asymptotic expansion 

^(t, r ) = V 3 > + ±F™ + ^(F- MF«) + 1 {^F^ - MF^j + 0(l/r A ) (20) 

where 0{l/r A ) means an expression whose absolute value is bounded by A/r 4 , with A 
some positive constant, F = F(t — r*), is an arbitrary real function, and F^(x) = 
d n F/dx n . We assume also that F(x), F^(x), and F^(x) are bounded, that F^(x), 
is square integrable, and that and 

lim x 2 F {3} (~x)F {4} (-x) = 
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(These are essentially a "finite energy", plus "smoothness at infinity" conditions). 

Replacing @, (|7|), and (||) in (|19D , and expanding the resulting expression after 
replacing the asymptotic expression (EH) for x/j, we find 



9* 



M 



J?{3} F {4} + (_p{3}\: 



+ 



5_ 

9* 



i^! i?{3} F {4} , ^fp{3}\2 , _2_ p{2} p{3} 

45 18 V ; 15 



+ 0(l/r) 



(21) 



From the assumptions made above on F, it is clear that the time integrals of the 
second, third and fourth terms on the right hand side of (|2i|) , required to obtain 
vanish in the limit of large r. We therefore have 



lim 

r— >oo 



dt 



dt 



1 

15^0 



(F^) 2 dt 



(22) 



which leads immediately to 
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